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Abstract 

We present a stochastic and variational aspect of the Lax-Friedrichs scheme applied 
to hyperbolic scalar conservation laws. This is a finite difference version of Flem- 
ing's results ('69) that the vanishing viscosity method is characterized by stochastic 
processes and calculus of variations. We convert the difference equation into that 
of the Hamilton-Jacobi type and introduce corresponding calculus of variations 
with random walks. The stability of the scheme is obtained through the calculus 
of variations. The convergence of approximation is derived from the law of large 
numbers in hyperbolic scaling limit of random walks. The main advantages due to 
our approach are the following: Our framework is basically pointwise convergence, 
not L 1 as usual, which yields uniform convergence except "small" neighborhoods of 
shocks; The convergence proof is verified for arbitrarily large time interval, which 
is hard to obtain in the case of flux functions of general types depending on both 
space and time; The approximation of characteristics curves is available as well 
as that of PDE-solutions, which is particularly important for applications of the 
Lax-Friedrichs scheme to the weak KAM theory. 
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1 Introduction 



The Lax-Friedrichs scheme is one of the oldest, simplest and most universal technique 
of computing PDEs. There is the huge literature on the stability and convergence of 
the Lax-Friedrichs scheme based on the L 1 -framework, particularly for shock waves. In 
this paper we investigate the Lax-Friedrichs scheme in terms of scaling limit of random 
walks and calculus of variations and present several useful new results, relating the 
scheme to principles of probability theories and calculus of variations. We refer to only 
the Lax-Friedrichs scheme applied to hyperbolic scalar conservation laws. However our 
results imply that finite difference methods applied to evolution equations are likely to 
possess stochastic or variational structures. It is a meaningful effort to investigate the 
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Lax-Friedrichs scheme or other schemes applied to PDEs of various types by similar 
approach to ours. 

We consider initial value problems of the inviscid hyperbolic scalar conservation law 



ut + H(x, t, c + u) x = inTx(0,T], 

u(x, 0) = u°(x) G L°°(T) on T, f. u°(x)dx = 0, 



where c is a parameter varying within an interval [cq, c\] and T := R/Z is the standard 
torus. The assumptions for the flux function H are the following (Al)-(A4): 

(Al) H(x,t,p) : T 2 x R^R, C 2 (A2) H pp > (A3) lim E^lllA = +QO , 

|p|->+oo \p\ 

By (Al)-(A3), we have the Legendre transform L(x, t, £) of i?(x, t, •), which is now given 
by 

L(x,t,£) = sup{£p- H(x,t,p)} 

pgr 

and satisfies 

(Al)' L(x, t,():T 2 xl^M, C 2 (A2)' % > (A3)' lim = +00 . 

|£|^+oo |4| 

The last assumption is 

(A4) There exists a > such that l-L^I < a(\L\ + 1). 

Throughout this paper, T-dependency is identified with IR-dependency with Z-periodicity 
and T with [0, 1). (Al) and (A2) are standard in the theories of conservation laws. (A3) 
is necessary, when we introduce a variational approach stated below to our problems. 
(A4) is used for derivation of boundedness of minimizers of the variational problems. 
The problems (11. lft appear not only in continuum mechanics but also in Hamiltonian 
and Lagrangian dynamics generated by H and L [6], pU], [5], [Hj. In the latter case the 
periodic setting is standard. We remark that the whole space setting is also available 
with additional assumptions for H required by variational techniques. 

It is sometimes very convenient to introduce initial value problems of Hamilton- Jacobi 
equations which are equivalent to (11.11) 



;i.2) 



v t + H(x,t,c + v x ) = h(c) inTx(0,T], 
v(x, 0) = v°(x) G Lip(T) on T, 



where h : [co, c\] — > R is a given function. We consider (11.11) and (II. 2p in the class of 
generalized solutions called entropy solutions and viscosity solutions respectively. Such 
solutions exist in C°((0, T]; L°°(T)) and Lip(T x (0,T]). If u° = v° x , then the entropy 
solution u of (II. ip and the viscosity solution v of (11.21) satisfy u = v x (u is considered as 
a representative element). From now on we always assume that u° = v®. 

One of the central achievements in the analysis of (11.11) and fll.2[) is that they are 
closely related to the deterministic calculus of variations: The value of v at each point 
(x, t) G T x (0, T] is given by 

(1.3) v(x,t)= inf ( rL c (7( S ), S ,7 , (5))^ + ^°(7(0))j + McK 

7 eAC, y(t)=x \J I 
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where AC is the family of absolutely continuous curves 7 : [0,t] — > M and L c (x,t,£) := 
L(x, t, £) — c£ is the Legendre transform of H(x,t,c + •) (see e.g. [2]). We can find a 
minimizing curve 7* of ( II. 3p . which is a backward characteristic curve of u,v and a C 2 - 
solution of the Euler-Lagrange equation associated with the Lagrangian L c . If the point 
(x, t) is a regular point of v (i.e. there exists v x (x, t)), then the value u(x, t) = v x (x, t) is 
given by 

(1.4) u(x,t)= t L c x ( 1 *(s),s n *'(s))ds + u°( 1 *(0)), 



Jo 

where u° is supposed to be rarefaction-free, or equivalently v° is semiconcave. (Otherwise 
m°(7*(0)) needs to be replaced with L£(7*(0), 0, 7*'(0)). See Lemma I3T21) We remark 
that, since v is Lipschitz, almost every point are regular and (II .4p is valid for almost 
every point. 

The representation formula ( 11. 3p is a strong tool not only in the analysis of fll.lj) and 
( II. 2p but also in many applications of them to other fields such as optimal controls [9] 
and dynamical systems [7]. It should be noted that the variational approach to ( fi.ij) 
and U.S\) based on U.3\) and ( |i.^| ) also contributes approximation theories of U.l\) and 
U. 2\) by the vanishing viscosity method and the finite difference method. The first case 
is announced by Fleming jH] and the latter case is the theme of this paper. 

First we recall the results of Fleming. Let us consider initial value problems of 

(1.5) + H(x,t,c + u u ) x = vu v xx , 

(1.6) < + H(x, t, c + <) = h(c) + vvl x {v > 0) 

with the same setting as (II. ip and (II. 2p . The solutions u u and v u are also related to 
calculus of variations which are not deterministic but stochastic: The value of v v at each 
point (x, t) G T x (0, T] is given by 



;i.7) v"(x, t) = inf E 



* L c (Y(s),s,C(Y(s),s))ds + v°(Ym 







h(c)t, 



where E stands for the expectation with respect to the Wiener measure and 7" is a 
solution of the stochastic ODE 



(1.8) dY(s) =C(Y(s),s)ds + V2udB(t-s), Y(t)=x. 

Here B is the standard Brownian motion. There exists the unique minimizing vector 
field of (11.71) . The value u u (x,t) is given by 



[1.9) u u (x,t)=E 



1 L c x (Y*(s), s, r(Y*(s), s))ds + uV (0)) 



where 7^* is a solution of ( 11.81) with £ v = . It is proved from a stochastic and 
variational point of view that, for v 0+, v u converges uniformly to v with the error 
0(y/u) and u v converges pointwise to u = v x except for points of discontinuity of u. 
In particular, u u converges uniformly to u without an arbitrarily small neighborhood 
of shocks. The proof indicates how the stochastic variational formula 1 11- 7ft and (II. 9ft 
tend to the deterministic ones (11. 3D and (II. 4ft . Asymptotics of fl 1.8ft for v — > plays a 
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central role. Fleming's approach yields much information and concrete pictures of the 
vanishing viscosity method. In particular we can see how the parabolicity disappears to 
be hyperbolic. 

The purpose of this paper is to establish a stochastic and variational approach to the 
finite difference method with the Lax-Friedrichs scheme. We discretize the equation of 
(11. ip by the Lax-Friedrichs scheme: 

q ~ ( m 2 m+2) , H(x m+2 , tk, c + m^ +2 ) — H(x m , tk, c + u k m ) _ 

V " / A J- O A „ 



We will see in the next section that we can find a difference equation which approxi- 
mates the equation of (jl.2p and is equivalent to (ll.lOp in the sense that = (v^ +1 — 
<-i)/2Ax: 

fc+1 _ ("m-l+^m+l) k _ k 

(1-11) ^ + H(x m ,t k ,c+ Vm+1 2/ ^ m - 1 ) = h(c). 

We present stochastic calculus of variations associated with (11.111) that we minimize the 
expectation of a discrete action functional among space-time inhomogeneous random 
walks in AxZ x AtZ. This yields representation formulas of and u k m similar to 
(jl.7p and (11.91) . The probability measures of such random walks are no longer related 
to the Winer measure. This is the main difficulty of our arguments: We need to study 
the asymptotics for Ax, At — > of the random walks generated by arbitrary transition 
probabilities depending on space and time, under hyperbolic scaling < A < At/ Ax < 
Ai. We will see that the continuous limit of minimizing random walks is deterministic. 
In other words, we obtain the law of large numbers, where the minimizing random walks 
converge to the minimizing curves for u, v. This proves convergence of the Lax-Friedrichs 
scheme. It is interesting to note that, under diffusive scaling Ax 2 /At = 2v > 0, the 
continuous limit of a certain class of random walks is the Brownian motion or diffusion 
processes, and the solutions of (jl.lOp and (11.111) converge to these of (ll.5p and (jl.6j) . 
Our approach also yields much information and concrete pictures of the finite difference 
method with the Lax-Friedrichs scheme. In particular we can see how the "parabolicity" 
due to numerical viscosity disappears to be hyperbolic in terms of the law of large numbers. 
Here we point out several advantages due to our approach: 

(1) Stability of the Lax-Friedrichs scheme and therefore convergence of the scheme up 
to arbitrary T > is derived from variational techniques. 

(2) The pointwise convergence of u k m to u — v x is proved. In particular this yields 
the uniform convergence except neighborhoods of shocks with arbitrarily small 
measure. 

(3) The uniform convergence of to v with an error 0{yAx) is proved from the 
stochastic viewpoint. 

(4) The approximation of (backward) characteristic curves of (II. ip and f ll.2[) and its 
convergence are verified. 

Approximation of entropy solutions with the Lax-Friedrichs (also with other schemes) 
is basically based on the /^-frameworks with a priori estimates, where Ax, At-independent 
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boundedness of both and its total variation must be verified. Our stochastic and vari- 
ational approach is totally different from the L 1 -frameworks and proofs are simpler. 

In the case of flux functions which are independent of x and t, many details are 
known. Crandall and Majda [I] prove stability and /^-convergence properties of mono- 
tone difference approximations in a rather general setting, where a flux function H(p) 
is not necessarily convex. Tadmor [17] shows the Lipschitz one-sided boundedness 
( u m+2 ~ u m)/2Ax < which guarantees time-global stability. §abac [13] proves 
that the optimal /^-convergence rate of u k m — > u is 0(V Ax). 

In the case of flux functions which depend on both x and t, the problem becomes 
much harder: Oleinik [12] extensively investigate the Lax-Freidrichs scheme in this case, 
where stability and /^-convergence are proved with restricted T > 0. This restriction is 
not satisfactory, because admissible T must be less than J^o^ loo yjpjdr w hL a monotone 
increasing C^-function V such that sup x t ^.j |p|< r \H x (x,t,p)\ < V(r) (V(r) = ar 2 is often 
the case, due to H pp > 0). Nishida and Soga [H] show the time-global stability of the Lax- 
Freidrichs scheme as well as the long time behavior that any difference solution converges 
exponentially to a periodic state as k — > oo yielding Z 2 -periodic entropy solutions. This 
argument is a generalization of results of Oleinik [12] and Tadmor [T7j. However they 
still assume that the flux function is of the form H(x,t,p) = |p 2 — F(x,t). It seems 
extremely hard to prove stability of the Lax-Friedrichs scheme for arbitrary T > with 
H(x,t,p) satisfying (Al)-(A4), in a similar approach. 

There are also lots of works on approximation of viscosity solutions of Hamilton- Jacobi 
equations. Many of them are done independently of entropy solutions. We remark that, 
even in the case of (jl. ip and (jl.2p . convergence results for viscosity solutions do not 
necessarily imply these for entropy solutions. We point out Crandall and Lions [3] and 
Souganidis [IB], where convergence results for viscosity solutions with general schemes 
are established as well as error estimates of 0(\/ Ax). Our results provide interpretation 
of the order V Ax from the stochastic viewpoint. 

Finally we refer to applications of our results to the weak KAM theory. The weak 
KAM theory makes clear the important connection among entropy solutions, viscosity 
solutions and the Hamiltonian dynamics generated by H(x,t,p). Approximation theo- 
ries of the weak KAM theory should be developed with approximation methods which 
provide all of entropy solutions, viscosity solutions and their characteristic curves at 
the same time, because the weak KAM theory requires three of them to connect the 
PDEs with important properties of the Hamiltonian dynamics. Fleming's results in- 
dicate that the vanishing viscosity method meets the requirement. Bessi [1] and other 
authors successfully exploits Fleming's approach to develop smooth approximation in the 
weak KAM theory. Nishida and Soga [IT] develop difference approximation in the weak 
KAM theory by the Lax-Friedrichs scheme with the usual /^-framework, where some 
arguments are not yet mathematically completed. Our results here indicate that the 
Lax-Friedrichs scheme also provides approximation to all of entropy solutions, viscosity 
solutions and their characteristic curves at the same time. The results of this paper can 
be strong tools for numerical analysis of the weak KAM theory. 
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2 Results 



2.1 Equivalent schemes for conservation laws and Hamilton- 
Jacobi equations 

Let N, K be natural numbers. The mesh size A = (Ax, At) is defined by Ax := (27V)" 1 
and At := (2K)' 1 . Set A := At/Ax, x m := mAx for m G Z and tk '■= kAt for 
k = 0, 1, 2, • ■ ■ . For x G M and t > 0, the notation m(x), /c(t) denote the integers m, A; 
for which a; G [x m ,x m + 2Ax),t G \tk,H + At). Let (AxZ) x (AtZ> ) be the set of all 
(x m ,t k ) and 

£w C (AxZ) x (AtZ> ), g odd c (AxZ) x (AtZ> ) 

be the set of all (x m , tk) with k = 0, 1, 2, • • ■ and m G Z with m + k =even, odd. We 
call Geven, Godd the even grid, odd grid. We consider the discretization of (11. ip by the 
Lax-Freidrichs scheme in Q even : 



„,Jfc+l ( u m+«m + 2 ) 



u 



m+1 2 



H(x m+2 , t k , c + u k m+2 ) - H(x m , t k , C + u k m 



(2-1) 



At 2Ax 



0. 



k — „.k 

m±2N u mi 



where for m = even 



(2.2) u° A (x) := — — / u°(y)dy for x G [x m - Ax,x m + Ax). 

Jx m -Ax 

Note that • 2Ax is conservative with respect to k and is zero for 

{m | 0<m<2N, m+k=even\ 

u° with the average zero. Now we consider a discrete version of (II .2D in ^o^: 



(2.3) 



fc + 1 _ ("m-l+^m+l) fe _ k 

^ + # (x m , t fc , c + m+1 2Ag m - 1 ) = 

W m+1 = v< h( x m+l)) V m+1±2N = V m+li 



where v\ is a function which converges to v° uniformly as A — > 0. We introduce the 
difference operators: 

fc+l _ VrHvfi ...k _ ...k 
D,w k+l — — 2 - D w k — m+1 m ~ 1 

In addition to the assumption u° = v®, we assume that 

(2.4) v° A (x):=v°(0)+ [ X u A (y)dy. 

Jo 

Note that u° A — > u° in L l and — > v° uniformly with || v A — v° u ° |U°° -3 Ax, as 

A — > 0. The two problems ( 12. ip and (12. 3p are equivalent under ( 12. 2 p and (12.4p : 

Proposition 2.1. Let u k m and v^ +1 be the solutions of \2.1\) and \2.$) with (CO)) and 
\2.1$ . Then the one is derived from the other. In particular we have D x Vm+i = u k m . 
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2.2 Random walks for the Lax-Friedrichs scheme 

We introduce space-time inhomogeneous random walks in Q 0( m, which corresponds to the 
stochastic ODE (11.81) in the vanishing viscosity method or characteristic curves of (11. ip . 

For each point (x n ,ti + i) G Godd, we consider backward random walks 7 which start 
from x n at ti + % and move by ±Ax in each backward time step: 

l = {l k }k=o,i,-,i+u l l+1 = x n , 7 fc+1 - 7 fc = ±Ax. 

More precisely, we set the following for each (x n ,^ + i) G Q m'- 

X k := {x m I (x m , t k ) G G odd , \x m - x n \ < (I + 1 - k)Ax}, k < I + 1 

g-.= |J (x k x {t k }) c g odd , 

i<k<i+i 

£:G3 (x m ,t k ) ^ e m e [-A" 1 , A" 1 ], A = At/ Ax, 
p : G 3 (x m , t k ) 1 y f m := - - -A£* G [0, 1], 

P ■ G 3 (x m ,t k ) H- p^ := - + -A^ G [0, 1], 

7 : {0, 1, 2, • • • , / + 1} 3 k ^ 7 fe G X fc , 7 ' +1 = x n , 7 fc+1 - 7 fc = ±Ax , 
Q : the family of 7. 

We regard p^ (respectively p^J as a transition probability from (x m , to (x m + Ax, t k — 
At) (from (x m , tk) to (x m — Ax, t k — At)). We control the transition of our random walks 
by £, which plays a "velocity field" -like role in G. We define the density of each path 
7 G Q as 

i<fc<i+i 

where p(7 fc ,7 fe_1 ) = p^*,-, (respectively p^ (7fc) ) if 7* - 7 fc_1 = -Ax (Ax). The density 
p(-) = pj(-; yields a probability measure of f2, namely 

prob(A) = ^2 ^(7; for ^Cfi. 

The expectation with respect to this probability measure is denoted by E^..^, namely 
for a random variable / : Q — > R 

^(■;C)[/(7)] :=X)M7;0/(7)- 

The above objects depend on the initial point (x n , U+i), but we omit to add the index to 
them for simpler notation. We remark that, since our transition probabilities are space- 
time inhomogeneous, the well-known law of large numbers and central limit theorem for 
random walks do not always hold in our case. Soga [15] investigate the asymptotics for 
A — > of the probability measure of Q under hyperbolic scaling, which will be important 
in this work. 
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2.3 Stochastic and Variational Representation of Approximate 
Solutions 



We represent the approximate solutions by a discrete action functional with the random 
walks and L c = L — c£. From now on we assume the following: 

Assumption. Suppose (Al)-(A4)- Let T > be arbitrarily fixed. The parameter c 
varies within [c ,Ci]. Initial datas are bounded: || u° — ||l°o< r, \\ v° ||c°< r - 

First of all we see a result, assuming also that there exists a solution u k m of fl2.ll) which 
satisfies the stability condition called the CFL-condition 

\H p (x mi t k ,c + u k J\ < X- 1 (A = At /Ax). 

Proposition 2.2. Suppose that we have the solution of 112.3]) for which := D x v^ +l 
satisfies the CFL-condition for all m and k = 0, 1, 2, • ■ ■ , k*. Then the solution is repre- 
sented for each n and 0<Z + l<fc* + l as 

(2.5) v l n +1 = ME^[ ^(T fe ,Vi,d (7 , ) )At + <(7 )]+/ i (c)t i+ i. 

0<k<l+l 

The minimizing velocity field £* exists, which is unique and given by 

C k + l = H p {x m) t k ,c + D x v k m+1 ). 

This proposition is informative, because a proof indicates how the Lax-Friedrichs scheme 
reveals the stochastic and variational structure. The proof also implies that Proposition 
12.21 holds only with the assumptions (A2) and (A3). 

Next we remove the assumption of the existence of v^+i satisfying the CFL-condition: 

Theorem 2.3. There exists X± > (depending on T, [cq, c±] and r, but independent of 
A ) for which we have the following: 

1. For any small A = (Ax, At) with A = At/ Ax < Xi, the expectation of the following 
functional with each n and < / + 1 < k(T) 

K +1 (0 E K , [ L c (7 fc ,^-i,e {7fc ))At + ^(7°)]+Mc)tz + i 

0<k<l+l 

has the infimum denoted by with respect to £ : G — > [— A -1 , A -1 ]. The infimum 
V^ 1 is attained by which satisfies |£*| < A^f 1 < A -1 . 

2. Define f^+i f or ea °h m and < k < k(T) as v^ +l := v%(x m+ i), := V^ +1 . 
Then the minimizing velocity field £* for each satisfies in G 

L^(x m , tfc, £, m ) = -D x V m+ i ^ £, m = H p (x m , tfe, c + D x v m+1 ). 

3. The above v k m+1 satisfies (TOP for0<k< k(T). 
This theorem immediately leads to one of our main results: 
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Theorem 2.4. There exists Ai > (depending on T, [co,Ci] and r, but independent of 
A) such that for any small A = (Ax, At) with A = At/ Ax < Ai we have the solution 
u k m of h2.1\) which is bounded and satisfies the CFL-condition up to k = k(T): 

\H p (x m ,t k ,c + u k m )\ < Af 1 < A -1 . 
Next we "represent" the solution of (12. ip : 

Theorem 2.5. For each n and < / + 1 < k(T), let £* be the minimizer for V^ +1 and 
7, //(•;£*) be the minimizing random walk for V^ 1 . Let be the minimizer for V^Xl an d 
7, //(•;£*) ^ e ^ e minimizing random walk for V^\. Then w^Vi satisfies 



v l+1 

u n+l 



v l+1 

u n+l 



< £ M , e) [ Yl ^(7 fc ,*fc-i,r™ ( y.))At + U i(7° + Ax)] + O(Ax) 

0<fc<Z+l 

^ %;£*)[ E ^(7 fc ,*fc-i,lt(7*))A* + «i(7°-Ax)] + O(Ax) 



0<fc</+l 

where O(Ax) stands for a number of (—9 Ax, 6 Ax) with 9 > independent of Ax. 



2.4 Convergence of Approximation 

We present convergence results of the stochastic and variational approach to the Lax- 
Friedrichs scheme. We always take the limit A = (Ax, At) — > under hyperbolic scaling 
< A < A = At/ Ax < Ai. We say that a point (x,t) 6Tx (0, T] is a regular point, 
if there exists v x (x,t). Note that regular points are nothing but points of continuity of 
u = v x and almost every points are regular. The minimizing curve of v(x, t) is unique, if 
(x,t) is regular (see e.g. [2]). 

Theorem 2.6. Let v be the viscosity solution of U.S\) and va be the linear interpolation 
of the solution v^ +1 of A2.3\) . Then 

V/\ — y v uniformly on T x [0, T] as A — > 0. 

In particular, we have an error estimate: There exists (3 > independent of A such that 

\\va-v ||c°(tx[o,t])< f^y/Ax. 

This result is consistent with the earlier literature cited in Introduction. The argument 
is based on the different viewpoint that the random walks become deterministic and our 
stochastic calculus of variations tend to the deterministic ones as A — > 0. 

Theorem 2.7. Let (x,t) E T x (0, T] be a regular point and 7* : [0,t] ->■ R be the 
minimizing curve for v(x, t) . Let (x n , ij+i) be a point of [x— 2 Ax, x+2Ax) x [t— At, t+At) 
and 7a : [0, t] — ^ M. 6e the linear interpolation of the random walk 7 generated by the 
minimizing velocity field £* for v 1 ^ 1 ■ Then 

7a — > 7* uniformly on [0,t] in probability as A — > 0. 

In particular, the average 0/7A converges uniformly to 7* as A — > 0. 
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The minimizing curve 7* is the genuine backward characteristic curves of u, v starting 
from (x, t). Therefore the Lax-Friedrichs scheme turns out to approximate not only PDE 
solutions but also their characteristic curves at the same time. 

Theorem 2.8. Let u = v x be the entropy solution of U.l\) and u& be the step function 
derived from the solution u k m of h2. namely ua(x, t) = u k m for (x, t) e [x m — Ax, x m + 
Ax) x [tk, t k + At). Then for each regular point (x, t) G T x [0, T] 

u^(x,i) — > u(x,t) as A — > 0. 

In particular, u& converges uniformly to u on (T x [0, T]) \ 0, where is a neighborhood 
of the set of points of singularity of u with an arbitrarily small measure. 

This convergence result is stronger than the one derived from the usual L 1 -framework in 
the following sense: The approximate solution ma converges pointwise a.e. and therefore 
"uniformly" to the particular representative element of the L x -valued entropy solution, 
which is the derivative of the corresponding viscosity solution and is represented as 



3 Proof of Results 

We prove the propositions and theorems in Section 2. 

Proof of Proposition ED It follows from (E2D and (E3} that D x v Q m+l = u° m . Sub- 
tracting D t v*+ X + H(x m , t k , D x v^ +1 ) = h{c) from + H (x m+2 , t k , D x v^ +3 ) = h{c) 
and dividing it by 2 Ax, we see that D x v^ +l is equal to u k m . 

Converse argument: We inductively define v^+i i n Godd as 

vt+i = vt-i + u k m ■ 2Ax, v k _ (m+1) = v k im _ 1} + u k _ m ■ (-2Ax) 

with v k := UqAx for k = even and Vq := for k = odd. Then v^ +1 satisfies D x v^ +l = u^, 
Vm+i±2N = ^m+i an d ^m+i = ^A^m+i) — u°(0). Using the equation of ( 12. ip and the 
equalities 

D t v k + X = D t v k+1 + {D t u k+l + D t u k+1 + ■■■ + D t u k ^\) 2 Ax for k = even, 
D t v k + l = D t v k+1 + {D t u k+1 + D t u k+1 + ■■■ + D t u k +\} 2 Ax for k = odd 

(note that D t VQ +1 = 0), we obtain 

D t v k m +1 + H(x m , t k , c + D x v k m+1 ) = P k , 

where for k = even, odd 

P k = D t v k+1 + H(x Q , t k , c + D x v k ), P k = D t v k+1 + H( Xl ,t k , c + D x v k ). 

We define v^ +1 as 

v° m+1 := v° m+ i + AO), v k m+1 := v k m+l - £ (P k ' - h(c))At. 

0<k'<k 

Then satisfies ( 12. 3p and D x v^ +1 = u k m . □ 
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Define fij^ := {7 6 O 7 fc = x m } and p k m := X! 7 en fe Mt)- We observe the following 
lemma, which follows from the definition of random walks. 

Lemma 3.1. 1. Pm(x) = 1- Hence {p k n ^} x &x k yields a probability of X k . 

xex k 

2- P k m = pfil'iSt 1 ! + ^i/t+i> ™here = z/x m±1 # X k +\ 

Proof of Proposition [2.21 Fix £ : G — > [— A^A^ 1 ] arbitrarily. It follows form the 
difference equation of (12.31) and the property of the Legendre transform that 



v l n +1 = \ V ' n+l - H(x n , t h c + D x v l n+1 )At + h(c)At 



= ■ (c + D x v l n+1 ) - H(x n , t h c + D x v l n+l )}At - c£ l n +1 At 

+{\ + \^ l n + >Li + (\- lK +1 K +1 + h(c)At 

< L c (x n , t h e 1 ) At + (i + hOLi + {\~ + He) At, 

where the equality holds, if and only if = H p (x n ,ti,c + D x v l n+1 ) G (— A _1 ,A -1 ' 
Similarly we have 

vL, < L^xn^t^ii^At + (1 + l -\e n ~iH;~2 + {\- l^LiX" 1 + h(c)At, 

v l n+ i < L c (x n+1 , e^ +1 ) At + (- + ^A^)^" 1 + (- - 2^+i) u n+2 + h(c)At, 

where the equality holds, if and only if £^ ±1 = H p (x n ±i, c+D x w^.\ +1 ) G (— A _1 ,A -1 ' 
Hence, by Lemma [3. 1[ we obtain 



< 



E ( E pt (X )L c (x, t k ^, e m[x) ) ) At + e PmiA 1 *) + ^x^+i - ^: 

k<l+l \x€X k J xex'- 1 



l<k<l+l 

Continuing this process, we obtain 



^< E ( E^)^^- 1 '^))) At+ E^" 



(x)^m(x) + M C )**+1- 

The equality holds, if and only if £^ = H p (x m , t^-i, c + ^i^+i) e ( — ^ _1 j We see 
that the first and second term of the right hand side, denoted by A\ and A 2 , are changed 
into 



A i = E { E ( E p(7;0)^(7 fc ,4-i,e (7fc) )}^ 

0<fc<« + l igj* 7 gQfe 

= E (E^(7;0^ c (7 fc ,t,-i,e (7fe )))At 

o<fc<m 7 en 

= £a*(7;0( E ^ c (7 fc ,^i,e (7fe) )At), 

7 €f2 0<fc<Z+l 

A * = E( E p(7;0)^( 7 o ) = E/ i (7^> 



^m( 7 )- 



£ is arbitrary and we conclude (I2.5p . □ 
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For simplicity, we introduce the following notation: £| 7 : k h-> £ 



771(7 



k ) ■ 



£(7,77) := / L c (7( S ), S ,r / ( S ))^ + ^ ( 7 (0)):AC([0,t])xL 1 ([0,t])^M, 
Jo 

UiliV) ■= Yl ^ C (7 fc ,4-i,^)At + ^( 7 ): J Dx J D^M, 

o<fc<m 

where D is the set of functions: {0, 1, 2, • • ■ , I + 1} — >■ M. 

Proof of Theorem 12.31 It follows from (A1)-(A4) that there exists ct\ for which we 
have \L C X \ < ai(|L c | + 1) for any c G [cq, Ci] and := | min{0, inf L c }\ is bounded. 

x,t,£,c 

Due to the lower boundedness of L c , has the infimum V^ +1 with respect to 

£ : G — >■ [—A" 1 , A^ 1 ] for any A. Note that, if £* is a minimizer for V^ +1 , then for each 
(%m,tk) E G V£ is attained by : G C G — > [— A" 1 , A^ 1 ] which is the restriction of £* 
to G whose vertex is (x m ,tk). Introduce 

a 2 := r + Tmaxi(x, t, 0), a 3 := cci{(l + 2L*)T + a 2 + r}, 
(3.1) Ai := maxjmax |if„(a;, t, Ci + 1 + r + a 3 )|, max t, c — 1 — r — a 3 )|}, 

9 :=T max L xx (x,t,£). 

Now we assume that, for each n and some I such that < / < k(T), the minimizer 
£* for V^ +1 satisfies |£*| < A^ 1 . This is true for / = 0, because of the boundedness 
l-Dx^m+il — r an d Proposition 12.21 We want to prove the following estimate under 
Ax9 < 1: 



(3.2) \D x V r 



x " n+2 1 



yl + l _ yl + l 

v n+2 v n 



2Ax 



< 1 + r + a 3 . 



Before proving, we observe the result from (13. 2\i : Consider £ l n ^_ 1 (^) with £ whose re 
striction to k < I + 1 coincides with the minimizers for V^ +1 and V^l Then, using; the 
property of the Legendre transform, we have for any 



+2- 

n+1 



(3-3) £^(0 = L^Xn^Mi^A^^^V^ + P^V^h^At 

> iC+i ■ (c + - c + D x V*Xl)}At - cCAAt 

+{\ + ^ l n \\)Vr + {\ - ^ffi)V££ + /i(c)At 

= " ^ n+2 - H(x n+1 ,t l+1 , c + DXXD At + He) At, 

where the equality holds, if and only if = H p (x n+ i, ti + ±, c + D X V^ 2 )- Therefore the 
minimizer £* for V^Xi exists and satisfies 

= H p (x n+l , t l+1 , c + D X V£1). 

In particular, we have |£*J[+x| < A]" 1 and |£*| < A7 1 because of (13.21) . Thus, by induction, 
we complete the proof of 1. and 2. 
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We prove (13.21) : Let £* be the minimizer for V^ +1 and /x, 7, fi, etc. be the notation of 
the random walks starting from (x n , Let £* be the minimizer for V^Xl and /i, 7, S7, 
etc. be the notation of the random walks starting from (x n+2 , We take £) with 

C( x m,tk) '■= C*( x m + 2Ax,tfc). Since is the infimum, we have 



£a(7,CI 7 )J + Hc)t l+1 , ViX\ = E K .. t) [£ A (7,rk)J + Hc)t l+1 . 
Since Q = {7 — 2Ax | 7 G f2} and /i(7 — 2 Ax; () = filj; £*), we have 

^n +1 < ^C;C)[^a(7,CI 7 )] +/i(c)t I+1 

= £>(7;?)[ E ^(7 fc -2Ax,t fc _ 1 ,ft ( - fe) )At + ^ ( - _ 2A J +/i(c)t i 



Therefore we obtain 



0</c<2+l 



(3.4) £>,t£g ^ 



£a(7,CI 7 ) 



2Ax 



*k > 
m(7' I )> 



0<fc<«+l 



-L c (7 fc -2Ax,t fe _ 1 ,r^))}At 



^(70) t; m(7 -2Aa!) 



2Ax 



0<fc<«+l 



O(Ax) 



Using (A4) and noting that £* is the minimizer, we see that the fourth line of (13.41) 
denoted by A is bounded from the below by — (1 + r + a 3 ): 



A > 



-%*)[ E «i(l + l^ c (7 fe ,^-i,et (W )l)Ai 



r — #A;r 



0<fc<Z+l 



0</c<2+l 



r — #A:r 



^ E «i(l + ^ c (7 fe ,4-i,f^)) + ^)At 



- axL^T - r - 1 



0<fc<Z+l 



«i% ; e)|^A(7,ri 7 ) +ai#A(-;e) w a(7°^ 



aiT — 2«iL ! „T — r — 1 



> -ai^(. ;0 ) >Ca(7,0|t 

> — «i«2 — air — QiT — 2a\L*T — r — 1 
= -(a 3 + r + l). 



air — a\T — 2a\L*T — r — 1 



Similar reasoning by £) with C(^m> ^fe) := CO^m — 2Ax, yields the upper bound 
D V l+1 

Since (13. 3p becomes an equality for the minimizing velocity field, we conclude 3. □ 



Theorem 12.41 follows from Proposition 12.11 and Theorem 12.31 Theorem 12.51 follows from 
(13.41) and the similar reasoning for the upper bound of D x ViX\. 
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We observe several properties of the function «:Tx [0, T] — > M defined as 
v(x,t):= inf £( 7 ,V) + McX v{-,0) :=v°(-), 

76AC, l(t)=x 

which is the viscosity solution of (jl.2p . There exists at least one minimizing curve for 
v(x,t), due to Tonelli's theory. 

Lemma 3.2. Let 7* : [0,t] — > R be a minimizing curve for v(x,t). 

1. The following regularity properties hold: 

L\tf{T),Trf'{T)) G a-^( 7 *(r),r) for < r < t, 
LI(y*{T),T,y*'{T)) G a>(7*(r),T) for < r < t, 

where d~v (d+v) is the subdifferential (superdifferential) . In particular v x (Y (r) , r) 
exists for <t <t and is equal to L|(7*(r), r, 7*'(r)). 

|7*'(t)| < A^f 1 /or < r < t, where \\ is given in \3.1\) . 

3. v is Lipschitz continuous. 

4- If {x, t) is a regular point, then we have for any < r < t 



u(x,t):=v x (x,t)= / L c x ( 1 *(s) 1 s ll *'(s))ds + LI( 1 *(T),T ll *'(r)). 
Jo 

If v° is semiconcave or u° is rarefaction-free, then L^(7*(0), 0, 7*'(0)) = u°( 7 *(0)). 
In particular, v is the viscosity solution of Al.fy) and u the entropy solution of M.l\) . 

Proof. 1., 3., 4. are known, but for reader's convenience we give a brief proof. 
1. It follows from the minimizing property of 7* that we have for any < r < t 

v(x } t) = J % L c ( 7 *(s), s, l*'(s))ds + u( 7 *(t), r) + h(c)(t - r) 

and 7*| [o,t] is a minimizer for 1^(7* (r),r). Hence for any e > such that r + e < t we 
have 

/r+e 
L c (7*(s), s, 7*'(s))ds + u(7*(t), r) + h(c)e. 

Define 7 : [0, r + e] ->■ R as 7(5) := 7*(s) + 1± p^ with 5^0 for s G [r,r + e] and 
7(5) := 7*(s) for s G [0, r], where 7* is a minimizing curve for f(7*(r) + 5, r). Then we 
have 

L c ( 7 (s), s, i(a))ds + v(j*(t) + 5, r) + /i(c)e. 

Therefore we obtain 

u( 7 *(r) + 5, t) - u( 7 *(r), r) > ^( 7 *(r), r, 7 *'(t))<5 + 0(e5) + O(y), 
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which means that for e = y/5 

, t,( 7 *(r) +5,t)- v(r(r), r) - Lg( 7 *(r), r l*\r))5 

limillt r-r. > 0. 

<5-)-0 \5\ 

Similar reasoning with 7 : [0, r] -> R defined as 7(5) := 7*(s) + s " (r " £) ^ for s G [r — e, t\ 
and 7(5) := 7 *(s) for s G [0, r — e] leads to the second inclusion. 

2. Since 7* is a minimizer, we have for 7 (s) = x 



L c ( 7 *(s), s, l*'(s))ds + u (7*(0)) < / ^ 5, 0)cfe + u (x) < a 2 . 

Jo 

Since 7* satisfies the Euler-Lagrange equation, we obtain 

m( 1 *(s),s, 1 *'(s))\ = |L!( 7 *(0),0,7*'(0))+ f Ll( 1 *(s),s, 1 *'(s))ds\ 



< r+ / a 1 (l + \L c ( 1 *(s),s ll *'(s))\)ds 

Jo 

< a 3 + r + l. 

3. follows from 2. 

4. We observe that for any < r < t and 5^0 

+ < J L c ( 7 *(*) + 5,s,7*'(s)) -L c ( 7 *(s), S ,7*'(s))^ 

+v(7*(r) + <y,r) -v(7*(r),r). 

Therefore, using 1., we obtain 



n(x,t) = ^.(x,t) < J L c x ( 1 *(s),s, 1 *'(s))ds + Ll( 1 *(r),r, 1 *'(r)), 

which holds for r = 0. 

For the converse inequality, we take a sequence Sj —> such that each (x + Sj,t) is a 
regular point of v. Let jj be the minimizer of v(x + 5j, t). Then we have lim ? -_ > . 00 v x (x + 
Sj,t) = v x (x, t). Therefore it holds that (7-,-, 7J) — > (7, 7') uniformly as j — > 00. A similar 
reasoning with 7 j yields 



M (x,t) = ^(x,t) > y ^(7*( S ), S ,7*'( S ))d S + L ? c (7*(r),T,7* , (r)), 

which holds for r = 0. We refer to [2] for the fact that v is the viscosity solution of (11.21) 
and u the entropy solution of (II. ip . □ 

We state a result of hyperbolic scaling limit of random walks. We introduce a random 
variable 77(7) = {rj k (7) }fc=o,i,2, -,1+1, 7 G f2 which is induced by a random walk with a 
velocity field £ and is defined as 

^+i := y+i ; ^fc( 7 ) := y+i _ ^ £(t fc ,, 7 fc ')At for <£;</. 

k<k'<l+l 
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Lemma 3.3. QlSl) Set ~o k := £ m( , [| 7 * - V k (l)\ 2 } and d k := E^y" - V k (l)\] for 
< k < I + 1. T/ien we have 

4.1+1 _ fk 

(d k ) 2 <~a k < Ax. 

A 

Remark. We always have a k — > as A — > for any £ under hyperbolic scaling. 
However the variance does not always tend to (see fTo^). This estimate is analogous 
to the following: Let rf be the stochastic process induced by the solution y u of M.8\) as 
rf'(s) = £"(s, 7 I/ (s)), rj u (t) = x. Then we have for s G [0, t] 

E[W{s)-Y{s)\} = V^E[\B{t-s)\}. 



Proof of Theorem 12.61 Since va and v is Lipschitz continuous, it is enough to show 
that \v l + x - v(x n ,ti +1 )\ = 0(y/Ax) for all < / + 1 < k(T) and n. Hereafter, ft, /3 2 , ■ • • 
are constants independent of Ax, At, n, I. Let 7* be a minimizer for v(x n , iz+i). By 2. of 
Lemma [3.21 we have |7*'| < A^ 1 < A" 1 . We take 

£(x m ,t k ) := 7*'(tfc), 

which yields a space-homogeneous random walk. Then we have 



v(x m t l+1 )=jr,{y*,y*') + h(c)t l 



+i> 



vi +1 < E 



Since £ is space-homogeneous, 77(7) and Ca(t){i)j CU) are independent of 7 and satisfy 

|7*(t fc ) -v k (l)\ < /3iAx forO<fc<Z + l, 
\£(Y,1*')-£a(v(i),Z\,)\ < /3 2 Ax+K;°(7*(0))-<(r / (7))|</3 3 Ax. 



By Lemma [3.31 we obtain 

v 1 ^ 1 -v(x n ,ti+i) < E^)\£ A (7,£\-,) - £a{v(i),€U) 



+ /3 3 Ax 



E K .#[ {L%l\t k ^ 1 *\t k ))-L\r ] \ 1 ),t k ^ 1 *\t k ))}At 



0<fc<Z+l 



+ (3 3 Ax 



< E.,1., 



£ /3 4 |7 fc -r/ fc (7)|At + r|7°-r / °(7)| 



0<fc<Z+l 



+ /3 3 Ax 



< /3 5 vAx. 



Let £* be the minimizer for v 1 ^ 1 , which yields a space-time inhomogeneous random walk. 
Let //a (7) : [0,£z+i] — > ffi be the linear interpolation of 77(7). Then we have 



,1+1 



E^*)^Ca(i,CU)\ +h(c)t l+ll 
v(x n , t i+1 ) < C{t]a{i), VAd)') + h(c)ti +1 for each 7 e Q, 
|jC(//a(7)> ^(7)') - (^(7)^*17) I < /^Ax for each 7 e ft. 
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By Lemma [3. 3[ we obtain 

v l n +1 -v(x n ,t l+1 ) > E K ,^\C A (j,CU)-C A (r]^),C\ 



-As Ax 



0<fc<Z+l 

+v A (l°)-v A (v°(l)) 



> -E*f)\ E Al7*-»? lb (7)|At + r|7 -»7 (7)l 



0<fc<«+l 



(3 6 Ax 



> -/3 7 VAx. 



□ 



In order to prove the convergence of minimizing random walks to minimizing curves, we 
observe the following lemmas: 

Lemma 3.4. Let 7* be the unique minimizer for v(x,t) . We define fore > 0, b > 
T s :={ 1 :[0,t]^R\ 1 eL i p, | T (t) - T *(t)| < e, | 7 '(»| < b, £(7, 7') < £(7*, 7*0 + 4- 
Then we have sup || 7 — 7* ||c°([o,t]) — ► 0, sup || 7' — 7*' ||z, 2 ([o,t]) — as e — > 0. 

Proof. Suppose that we have a sequence £j — > for which there exists 5 > such 
that sup 7gr Ej || 7 — 7* ||c°> 25 for all j. We can take a sequence 7, G T £j for which 
|| 7j — 7* ||c°> 8 holds for all j. We have a sub-sequence of jj, still denoted by 7^, such 
that 7j — >■ 7 uniformly. 7 satisfies || 7 — 7* ||c°> ^, 7(0 — 7*(0 — x anci £(7)70 = 
£(7*, 7*'). Therefore 7 is another minimizer for v(x,t), which is a contradiction. 

Set 5(e) := 7*(t) — 7 (t) for 7 G T e . Since 7* is the minimizer, we already have 
£(7*,7*0 < C(j + 6(e) n ') and \£(rf + 6(e), f) -£(7, V)l < fte (A > 1) for any 7 G P. 
Hence we have £ > £(7,7') — £(7*, 7*') > —j3%e. Therefore we obtain 

> |£( T ,7 / )-£(7*,7* / )l 

{L c (l(s), s, 7'(s)) — L c (7*(s), s, 7'(s)) 

+L C ( 7 *( S ), s, y(a)) - L c (7*(*), s, 7*'W)}da + «°(7(0)) - t>V(0)) 

> 

-^9 II 7 - 7* He" 
= ^(7*( S ), S ,7^))(7(^)-7*( 5 ))IS- 



' {^( 7 *(s), S , y'WXYto - 7*'(*)) + ^« " (Y00 " 7*'(*)) 2 }^ 



f^ l -L^.{ 1 '{s)- 1 *'{s)fds 



Pb || 7- 7* lie , 



where we remark that L^(7*(s), s, 7*'(s))' = L£(7*(s), s, 7 *'(s)) due to the Euler-Lagrange 
equation. Since is bounded away from zero on each compact set, we see that for con- 
stants /?io > 0, > 



8 s > /3io||7 / -7* / |ll 2 -/3iill7-7*llco. 
We have already proved that || 7 — 7* \\c°— > as e — > 0. This finishes the proof. 



□ 
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Lemma 3.5. Let f : [0, t] — > R be a Lipschitz function with a Lipschitz constant 9 
satisfying fit) = 0. Then we have || / ||co([o,t])< II / IU 2 ([o,t]) +\/\\ f IU 2 ([o,t])- 



Proof. We suppose that || / ||c*°> a/II / \\l 2 - Otherwise we get the conclusion. For 
I :={se [0,t] | > y/fThfi, we have || / ||| 2 > / 7 |/(s)| 2 rf S >|| / \\» \I\ and 

\I\ <|| / ||l 2 - The set I necessarily contains such that |/(s*)| =|| / \\c°- Since fit) = 0, 
there always exists s > such that f(s ) = ^\\ f \\ L 2 and [s*, s ] C /. Hence we have 
|/(s ) - /(s*)| < #|s - s*| < 0|I|. Therefore we obtain |/(s*)| < 0|/| + |/(s )|- □ 

Proof of Theorem 12.71 Consider ft A := {7 G H| || 7a — 7* [jco< e}, where 7a is the 
linear interpolation of 7 extended to s G [0,t]. We show that, for each arbitrarily fixed 
e > 0, we have profe(ft A ) — > 1 as A — > 0. 

First of all we prove the following: Introduce 



1 a •- 



{7 e ft I II ??a(7) - 7* ||c°< £ and || ^(7)' ~ 7*' IU 2 < e}, 



where ?7a(7) is the linear interpolation of 77(7) extended to s 6 [0, t]. Then, for each 
arbitrarily fixed e > 0, we have prob(T A ) — > 1 as A — > 0. In fact, we observe that 



.1+1 



E, 



^A(»7(T),ri 7 )J +h(c)t l+1 + 0(VA^) 
C(va(i),Va(i)')} + /i(c)t m + 0(v / A^), 



v(x, t) 



E„.u 



^a(7)^a(7)0 ->C(7*,7*')] +0(VA^) 
= 0(VAx). 
Since 7* is the minimizer, we have for all 7 G ft 

< £(t?a(7) + x-x n , 7^(7)') - £(7*, 7*') = £(^(7), ^(7)') - £(7*, 7*0 + P12 Ax. 
Consider ft+ := {7 G ft | >C(?7a(t),»7a(t)') - £(7*, 7*') > Axi}. Then we obtain 

£(7 7a (7)^a(7) , )-A7*,7* , )1 > E M7;r)Axi-/3 12 A 3 



O(VAx) 



76^ 



Therefore prob{Vt + ) = 0(Ax*). It follows from Lemma [3.41 that for Ax <C e 

ft \ ft + c r A . 

Thus we conclude that 1 — prob(Q + ) < prob(T e A ) < 1 and prob(T 6 A ) — > 1 as A — > 0. 

Since 77(7) converges to 7* uniformly in probability, for any e' > there exists 5 > 
such that if |A| < 5 we have 

^(■;D[l7 fc - 7*ftb)ll < ^(-;e)[l7 fc - + E»U*)l\v k (l) ~ 7*ftb)|] < e' for all fc, 

^ ( ,e)[||7A-7* M < e'. 

Define ft++ := {7 G ft | || 7a - 7* IU 2 > Ve>}- Since e' > £ M , e )[|| 7a - 7* IIl 2 ] > 
prob(£l ++ )\/~P, we have prob(£l ++ ) < ye 1 . Take e' = 0(e 4 ). By Lemma 13.51 we have 
I 7a — 7* ||c°< £ for all 7 G ft \ ft ++ . Therefore we conclude that ft \ ft ++ C ft A , 
1 - \fe' < prob(tt £ A ) < 1 and prob(tt £ A ) 1 as A ->■ 0. □ 
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Proof of Theorem [2781 Let (x n ,ti + i) be such that x G [x n — Ax,x n + Ax), t G 
[ti + i,ti +2 ). It is enough to show that u r ^} 1 — > u(x,t) as A — > 0. We set 



^(75 7') 



^x(7(s),s,7'(s))c/s. 



First of all we assume that v° is semiconcave, or u° is rarefaction- free. Let £* be the 
minimizer for v„. It follows from Theorem 12.51 that 



< [^'(^(7); ^a(7)0 + ^a(7° + Ax)J + O(VA^), 

where 77a(t) is the linear interpolation of 77(7) extended to s G [0, t]. Therefore, by 4. of 
Lemma 13.21 we obtain 

u l + + \-u(x,t) < JV^VM^aMO-^V^*') 

+u A ( 7 ° + Ax) - u°(7*(0))] + O(VAsc), 

1 



<( 7 u + Ax)-n u ( 7 *(0)) 



2Ax 



(7°+Az)+Ax 



«°(l/)-« (7*(0))^. 



(7°+Ax)-Ao; 



Since || r/ A (7) - 7* ||cr — > 0, || r} A (j)' - 7*' \\l 2 ^ 0, || 7a - 7* ||c?o— >- as A ^ in 
probability and w°(y) — >■ u°(7*(0)) as y — >• 7*(0), we obtain 

(3.5) lim sup(u l r ^ 1 — u(x,t)) = 0. 

A^O 

Similar reasoning with the converse inequality in Theorem 12.51 yields 

(3.6) hminf(< + + 1 1 - M (x,t)) = 0. 

We remove the assumption that v° is semiconcave, or u° is rarefaction-free. Hence we 
do not always have w°(7*(0)) (in this case, u° jumps up at x = 7*(0)). By 4. of Lemma 
13.21 and Theorem 12. 5[ we see that for any < r < min{tj + i,t} 

' L c (7*(s), s, i»(a))d8 + ^( 7 *(r), r, T *'(t)), 



u(x, t) 



m(7 fc ( T ))+l 



O(Ax) 



k(r)<k<l+l 



E 



m(7 fc (' r )) / 



O(Ax) 



fc(-r)<fc<«+l 



= E rt-,t*)[j L c x (r] A (-f)(s),s,TiA^)\s))ds 

+^(^(7)(r),r,r /A ( 7 ) , (r))] +0(v / A^). 

The average of ?7a(7)', denoted by 77^, converges to 7*' in L 2 . Therefore for any e > 
there exists 5 > such that if |A| < 5 we have || r/ A — 7*' IU 2 < e - I n particular we have 
< r* < min{ti + i,t} such that |77a( t *) ~7*'( r *)l < e - Taking r = r*, we conclude (13. 5p . 
Similar reasoning yields (13. 6p . □ 
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